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SYMBOLS 


r,  e 

u,  V 
P 

4^ 

R,  P.  a 

p.,  V 

P 

T 

6 

1  i! 

M 

Re 

Pr 

Pe 

\ 

T1 

h 

V 

S 

k 

U3 

q 


polar  coordinates 

radial,  tangential  velocity  components 

pressure 

streamf  unction 

see  Fig.  1 

absolute,  kinematic  viscosity 

mass  density 

temperature 

boundary  layer  thickness;  also,  "variation  of" 

"jump"  in  quantity;  "+"  refers  to  free  shear  layer 

Mach  number 

Reynolds  number 

Prandtl  number 

Pe'clet  number 

ratio  of  specific  heats 

/'P  ^ 

tjj  neck 
static  enthalpy 
velocity  vector 
specific  heat 
thermal  conductivity 
vorticity 
heat  flux  rate 


Subscripts 

r,  e 

partial  derivatives 

1.2,  3,  4, 

flow  field  subregions  (Fig.  1) 

DSL 

dividing- streamline  quantity 

D 

dividing  streamline  quantity  at 

t,  o 

stagnation  condition 

w 

wall  condition 

« 

dimensional  quantity 

H 

inviscid  flow  conditions  at  v  = 

B 

recirculation  region  quantity 

n 

component  normal  to  shock 

Superscripts 

f 

ordinary  derivative 

1 

first  approximation 

non-dimensional  quantity 
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ABSTRACT 


An  analytical  model  of  the  two-dimensional  laminar  near  wake  is 
described.  All  relevant  physical  features  of  the  flow,  including  the  boundary 
layer  expansion  at  the  shoulder,  free  shear  layers,  recirculation  regions 
and  recompression  region  are  included.  A  tractable  problem  is  formulated 
by  matching  approximate  solutions  for  each  of  these  regions  along  mutual 
boundaries.  A  set  of  coupled  algebraic  equations  is  derived,  and  numerical 
results  obtained  for  the  conditions  of  a  wind  tunnel  wedge  experiment.  Satis¬ 
factory  agreement  is  obtained  between  the  measured  and  theoretical  variation 
of  base  pressure  with  Reynolds  number.  Additional  computations  are 
carried  out  for  a  series  of  wedges  of  different  length,  apex  angle  and  wall 
temperature.  The  variation  of  base  pressure,  wake  angle  and  neck  enthalpy 
ratio  with  altitude  is  obtained  for  these  bodies.  It  is  shown  that  the  near 
wake  dimensions  scale  approximately  with  body  size,  and  that  the  neck  en¬ 
thalpy  ratio  has  a  significant  variation  with  body  size,  wall  temperature  and 
altitude.  The  effect  of  free -stream  velocity  on  the  neck  enthalpy  ratio  is 
seen  to  be  relatively  unimportant,  however.  A  number  of  extensions  of 
the  first-order  model  are  suggested. 
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Introduction 


The  laminar  near  wake  of  a  slender,  sharp-edged  hypersonic  body 
is  typical  of  a  large  class  of  unsolved  problems  in  fluid  mechanics  known  as 
separated  flows.  Typically  a  region  of  large  pressure  gradients,  reversed 
flow,  large  variations  in  Mach  number  and  non-adiabatic  conditions,  a  sep¬ 
arated  flow  probably  requires  the  solution  of  the  complete  Navier -Stoke s 
equations  for  its  accurate  description.  For  most  boundary  conditions, 
however,  this  problem  becomes  intractable.  Analytic  solutions  of  the  near 
wake  must  therefore  be  based  on  simplified  models,  which  then  yield  only 

approximate  results.  The  importance  of  this  problem^  and  the  complexity 

2 

of  obtaining  numerical  solutions  nevertheless  justify  such  an  approach. 

The  usefulness  of  a  flow  model  is  a  function  of  both  the  physical  content  of 
its  assumptions  and  its  capacity  for  refinement.  One  such  model,  con¬ 
sisting  of  several  subregions  of  the  flow  field  matched  along  common  bound¬ 
aries,  will  be  examined. 

The  flow  field  is  represented  in  Fig.  1.  It  is  described  in  terms 

of  inviscid,  boundary -layer  (free  shear  layer)  and  recirculation  regions. 

When  the  boundary  layer  thickness  is  very  small  compared  to  the  base 

height,  i.e.,  at  high  laminar  Reynolds  numbers.  Chapman's  original  theory^ 
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as  modified  by  Denison  and  Baum  should  be  adequate  for  two-dimensional 

5 

flow.  At  the  other  limit,  the  wake  of  a  flat  plate  or  a  needle,  Goldstein's 
and  Viviand  and  Berger's^  theories  neglect  the  "Base  Region"  completely. 
The  present  analysis  will  attempt  to  deal  with  conditirns  for  which  the 
boundary  layer  is  a  significant  fraction  of  the  base  height,  and  where  coup¬ 
ling  of  the  free  shear  layers  and  recirculation  regions  is  important.  In¬ 
clusion  of  the  recirculation  regions,  therefore,  constitutes  the  essential 
improvement  on  previous  theories. 

The  "inviscid"  regions  are  assumed  to  be  described  by  the  Prandtl- 


Meyer  relation,  and  the  free  shear  layers  (in  general)  by  the  Integral 
Boundary  Layer  Equations.  The  solution  of  the  free  shear  layer  with  this 

7 

method  by  Kubota  and  Dewey  suggests  the  approximation  of  a  straight-line 
velocity  profile  between  the  dividing  streamline  and  inviscid  flov  after  a 
short  transition  distance  from  the  shoulder.  The  relatively  hot,  low -density 
low -velocity  gas  occupying  the  recirculation  regions  suggests  that  they  can 
be  approximately  described  with  a  streamfunction  expansion  in  Reynolds 
number,  the  first  term  of  which  satisfies  the  biharmonic  equation.  An 
"elliptic"  system  which  letains  the  highest  order  derivatives  of  the  Navier- 
Stokes  equations,  the  biharmonic  equation  is  both  physically  and  mathe¬ 
matically  appropriate  in  these  regions.  The  influence  of  the  boundary  con¬ 
ditions  in  this  region  of  reversed  flow  naturally  results  in  the  formulation 
of  a  boundary  value  problem,  and  the  patching  requirements  along  dividing 
streamlines  demand  retention  of  highest-order  derivatives.  This  system 
allows  for  pressure  gradients  in  all  directions  and  can  be  solved  indepen¬ 
dently  of  the  energy  equation.  Most  important  of  all,  it  is  linear.  While  a 
streamfunction  obtained  from  the  biharmonic  equation  is  a  valid  represent¬ 
ation  of  a  low  Reynolds  number,  bounded  flow  field,  an  examination  of  the 
actual  limitations  of  this  approximation  is  given  in  the  next  section. 

The  requirements  for  matching  the  subregions  are  the  continuity  of 
pressure,  velocity  and  shear.  To  make  the  problem  tractable,  some  as¬ 
sumptions  about  the  shape  of  the  boundaries  of  each  region  are  made.  It  is 
assumed  that  1)  the  flow  is  symmetric  about  the  center-line;  2)  the  dividing 
streamlines  are  straight  and  3)  the  displacement  thickness  of  the  free  shear 
layers  is  constant.  Since  the  change  in  free  shear  layer  thickness,  6^,  with 
distance  is  relatively  small  when  the  expanded  boundary  layer  thickness, 

6'2f  is  a  significant  fraction  of  the  base  height,  the  assumption  that  6^^ 
is  justified.  Thus,  we  are  assuming  a  constant  pressure  inviscid  region 
(until  the  flow  turns  near  the  center-line)  and  cannot  hope  to  match  pressure 
at  more  than  two  points:  at  the  base  wall  and  at  the  important  rear  stagna¬ 
tion  point.  Furthermore,  by  avoiding  the  solution  of  the  free  shear  layer 
we  are  restricted  to  matching  of  velocity  and  shear  at  one  point  on  the 

*  M.  Bloom  has  pointed  out  to  the  author  that  Ting  and  Bloom®  first  sug¬ 
gested  this  approach  for  cavity  flows. 
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dividing  streamline.  The  near  wake  geometry  is  then  completely  specified 
by  the  free  shear  layer  thickness  and  the  wake  angle,  p.  The  matching  con¬ 
dition  on  the  dividing  streamline  is  thereby  reduced  to  finding  the  magnitude 
of  the  dividing  streamline  velocity  which  is  consistent  with  continuity  of 
shear.  Finally,  the  wake  angle  is  determined  by  the  continuity  of  pressure 
at  the  stagnation  point  of  the  dividing  streamline,  the  well-known  "Chapman 
recompression  condition".  Calculation  of  the  stagnation  pressure  requires 
knowledge  of  the  temperature  field,  and  this  is  determined  by  an  approxi¬ 
mate  solution  of  the  energy  equation.  From  this  solution,  it  is  possible  to 
calculate  the  stagnation  enthalpy  of  the  dividing  streamline  as  a  function  of 
all  the  flight  parameters.  Since  the  dividing  streamline  stagnation  point  is 
likely  to  be  the  hottest  point  in  the  flow  field  of  a  sharp,  slender  body,  this 
calculation  is  of  considerable  practical  importance. 

I.  Approximate  Solution  of  the  Recirculation  Region 

The  "vorticity  diffusion  equation",  the  reduced  form  of  the  incom¬ 
pressible  Navier-Stokes  equations,  is  written  in  non-dimensional  form  as 

-T.  V(V^4/)  =  0  (1) 

let  us  assume  il/  to  be  an  analytic  function  of  lie  and,  therefore,  an  expan¬ 
sion  of  the  form 

=  4*q  Re4i|  +  Re^4'2  +  .  .  .  • 

which  leads  to  a  velocity  field 

V  =  V  +  ReV ,  +  .  .  .  . 
o  1 

The  above  expansion  satisfies  the  equation  of  motion  for  arbitrary  Re  only 
if  coefficients  of  all  powers  of  Re  vanish  simultaneously:  these  conditions 
are 

=  0 

II  If 

etc.  etc. 
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where 


The  expansion  proposed  is  expected  to  be  valid  up  to  some 

an  upper  bound  to  Re  is  defined  by 

max  ^ 


The  first  term  in  the  expansion  satisfies  the  biharmonic  equation  (the  Stokes 
equations),  and  when  substituted  back  into  the  full  momentum  equations,  the 
first  order  pressure  distribution  may  be  calculated.  For  example,  along 
the  dividing  streamline  in  a  Cartesian  coordinate  system  constructed  as 
shown  in  Fig.  1, 


1 


u 


u 

o  x 


o 


and  the  important  inertia  terms  appear.  In  fact,  available  experimental 
17 

evidence,  together  with  the  solution  that  will  be  obtained  here  (Fig.  3), 
indicate  that 


1 

Ri 


over  most  of  the  recompression  portion  of  the  dividing  streamline  (velocity 
profile  inflection  point  is  at  =  0)  for  u^R/v^  >  10.  During  recompression, 
then, 

P.,.  ~  P*  constant.  , 

OX  t 

DSL  °  DSL  ° 


The  assumption  of  a  constant  total  pressure  recompression  process  will 
therefore  be  made.  It  is,  of  course,  possible  to  carry  out  the  integration 
to  obtain  _ 


P  —  P 
stag  base 
®o  o 


/^r '  1 

/  “o  -  "o-x 

•4  L  °J] 


and  obtain  a  higher  stagnation  pressure.  This  effect  is  obviously  Reynolds 
number -dependent. 

It  is  important  to  note  that  the  usual  test  of  the  validity  of  the 
Stokes  Streamf unction",  namely  that  the  ratio  of  inertia  terms  to  viscous 
terrhs  be  required  to  be  less  than  unity,  has  not  been  used.  Weinbaum^ 
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has  pointed  out  that  this  commonly  used  criterion  is  too  severe.  Wherever 

this  test  leads  to  ratios  less  than  unity,  however,  the  expansion  used  above 

will  be  valid.  Since  the  calculation  is  simple,  we  will  compute  this  ratio  in 

the  radial  and  tangential  momentum  equations; 

^  2 
ruu  +  vu,  —  V 
T  6 

— - 1 — 71 - 

rp  +  “r-?"-^F"ee- 


S  .  i 

O  V 

These  may  be  written  (with  r  =  r/R,  u  =  u/uj^) 
S 

-fe-  =  G^(r,0.p) 

S 

-fe-  = 


.V  .  uv 
uv  +  —  V  •  + - 

T  r  0  T 


T - 

—  V 

r  r 


"T - —2 - 

^2  '"ee 


rr 


r 


The  equations  of  motion  are  written  in  polar  coordinates  whose  origin  is  at 
the  stagnation  point  (Fig.  1),  the  resulting  concavity  in  the  base  wall  being 
insignficant  for  small  wake  angles.  The  streamfunction  is  then  defined  by 


V  =  lb 

r 


(2) 


The  first  term  in  the  expansion  in  Re,  must  therefore  satisfy 


l^'^rrrr  ^  r  ^rrr  ^‘^rr  ^  ^ ‘^^r  '^9699 

+  7^00rr  -  4  '<'00r  +  ^'^99  ]  "  ® 


(3) 


The  boundary  conditions  are: 
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(i) 

at  r  = 

0  : 

u  = 

0  =  V  ' 

(ii) 

on  r  = 

R  : 

u  = 

0  =  V  1 

(41 

(iii) 

on  0  = 

+p: 

u  = 

V  =  0  1 

1 

Note  that  must  be  determined  from  the  matching  conditions,  and 

the  boundary  conditions  reduce  to: 


(i) 

at 

r 

=  0  ; 

?+«  =  »  = 

!  (5) 

(ii) 

on 

r 

=  R  ; 

+^  =  0  =  4,^ 

(iii) 

on 

0 

=  _+p: 

-F+e  =  “dSl'*''’  ‘t'r  '  ° 

It  is  known  that  the  solution  to  such  a  Boundary  Value  Problem  exists  and  is 
unique.  It  is  further  known  that  the  solution  satisfies  (in  addition  to  the 
Boundary  Conditions)  the  Variational  Principle: 

V^^ip^eilj^rdrde  •=  0  (6) 

The  Ritz  Method  may  be  used  to  obtain  an  approximate  solution  for  4'^(r,0). 
Expanding  4*^^  in  a  set  of  suitable  functions  which  satisfy  the  Boundary 
Conditions; 


^o 


Ea  g  (r)f  (0) 
,  m®m  m 
m=  1 


(7) 


Ifboth  g  and  f  are  assumed,  based  on  the  expected  behavior  of  the  solu- 
^m  .m  ^ 

tion,  substitution  into  Eq.  (6)  will  result  in  an  infinite  set  of  algebraic  equa¬ 
tions  for  the  a  .  A  more  accurate  procedure  is  to  assume  only  the  e  Ir)- 
functions,  and  solve  ordinary  differential  equations  for  (which  include 

the  3-^).  Accurate  representation  of  the  velocity  profile  will  be  important 
in  the  matching  procedure,  and  this  latter  approach  entails  little  difficulty. 
Substituting  Eq.  (7)  into  (6),  and  noting  that  ^f  is  arbitrary,  yields; 


m=  1  •'o 


g  V^(g  f  )dr  = 
®p  ®’m  m 


(8) 


for  p  =  1,  2,  3,  ...  . 
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an  infinite  set  of  coupled,  ordinary  differential  equations  with  constant  co¬ 
efficients.  The  general  solutions  are 

f  =  c  ,f  ,  +  c  _f  _  +  c  -f  -  +  c  .f  .  (9) 

m  ml  ml  m2  m2  m3  m3  m4  m4 


with  the  c  's  determined  by  the  remaining  boundary  conditions  on  0  =  +  p. 
If  the  are  chosen  to  be  an  orthonormal  set,  important  simplifications 

are  possible.  In  addition,  we  may  then  expand  such  a  set  and 

obtain: 


m=  1 


where  b 


(If  non -orthogonal  functions  are  used,  Eqs.  (8)  are  coupled;  in  addition,  the 
b^  must  be  determined  by  another  procedure). 

The  symmetry  of  the  flow  may  now  be  used  to  eliminate  the  even 
functions  in  f^(since  v|i  must  be  odd;  i.  e.,  anti -symmetrical)  and  obtain 

^m  “  *^ml^ml  ^  ^m2^m2 
The  conditions  [(5)-(iii)]  may  be  written; 


g  f  (P)  =  0 

®m  m 


Substituting  Eq.  (11)  into  (12)  yields 

for  m  =  1,  2,  3,  .  .  .  . 


♦calculated  for  various  geometries  by  Chandrasekhar;^  "completeness' 
has  not  been  proved  for  this  or  any  other  applicable  set. 
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This  system  must  be  truncated  at  some  m  =  M,  giving  a  total  of  2M  equa¬ 
tions  for  2M  constants.  But  UQgj^(r)  can  then  be  matched  at  only  M  points 

to  determine  the  M  values  of  b  .  It  should  be  noted  here  that  no  proof  can 

m 

be  given  for  completeness  of  the  g^  set,  and  the  convergence  of  as  M-*oo 
can  only  be  shown  by  "numerical  experiment". 

For  these  boundary  conditions,  orthonormal  functions  have  the  dis¬ 
advantage  of  being  difficult  to  integrate,  and  it  is  helpful  to  choose  another 
"relatively  complete  set"  suggested  by  Kantorovich, 


g  =  r 
®m 


l-il 


In  addition  to  satisfying  the  boundary  conditions  [(5)-(i)  and  (5)-{ii)],  these 
functions  have  the  property  that  Eqs.  (1)  are  bounded  as  r  —  0, 

Choosing  m  =  p  =  M=  lin  Eq.  (8)  yields  an  equation  for  f  jt 

fj""  +  15.  33f^"  +  75.66f^  =  0  (15) 

In  this  first-order  approximation,  we  obtain 


2  fi  ,  , 

=  -r  1-5^  f,' 


Defining 


“d  ■  “dsl 


we  calculate  from  Eq.  (16) 


“DSL  = 


(^)  [^'k] 


The  conditions  to  be  satisfied  by  Eq.  (15)  are  then 

fl(P)  =  0  J 


and  the  solution  is 
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(r.fl)  = 


,3„  „2 


r  (1— r)^(z  sinh  a.0  cos  hQ  —  cosh  SlO  sin  b0) 


where 


r  =  b  cosh  a(3  cos  bp  +  a  sinh  ap  sin  bp 

—  z  (a  cosh  ap  cos  bp  —  b  sinh  ap  sin  bp) 


tan  b£ 


tahh  ap 


;  a  =  .  74,  b  =  2.  86 


and  sJ'^^/Re  may  now  be  evaluated  for  representative  wake  geo- 
metrics.  In  particular,  calculations  carried  out  at  121  mesh  points  for 
p  =  .  10,  .  15,  .  20,  .  30,  .40  and  .  50  radians  yield  the  following  results: 

a)  the  variation  of  S  and  S  with  P  is  small 

o  o 

»>  .2 

“Re - "Re“  general 

(r) 

c)  /Re  varies  many  orders  of  magnitude  over  the  recircu¬ 

lation  region,  but  a  "histogram"  display  indicates  the  over¬ 
whelming  (80%)  frequency  of 

S 

^  (y(5x  10“^) 

Thus,  this  severe  test  results  in  being  a  good  representation  of  ip  over 
80%  of  the  region  for  Re  =  p^Uj^R/jU^  up  to  200.  Since  it  can  be  shown  that 
typical  values  of  Re  are  of  this  size,  the  expansion  is  valid  over  most  of  the 
region.  Finally,  since  ip^  satisfies  the  boundary  conditions  exactly,  vp^  is 
an  exact  solution  on  the  dividing  streamline,  precisely  where  the  more 
severe  comparison  of  various  derivatives  of  ip  would  indicate  ip^  to  be  a  poor 
representation  of  Thus,  ip^  should  be  a  satisfactory  engineering 

estimate  of  vp  for  the  low  Re  cases  under  consideration.  Furthermore,  two 
additional  refinements  can  be  made: 

A.  Calculation  of  di, 

^  12 
Using  the  variational  principle 


5 -  2vp^ 


Vo*  V(V  ip^)>  jrdrd0  =  0 
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we  can  expand  4jj 


as 


oo 

m=  1 

and  proceed  as  before  for  imposing  homogeneous  bound¬ 
ary  conditions  on  f  ,  and  g 

ml  *ml 

B.  Approximation  of  the  Inertia  Terms  with  a  Modified  Oseen 
Technique; 

Since  the  high  Reynolds  number  recirculation  region  has  been 

23 

shown  by  Batchelor  to  consist  of  an  "inviscid  core"  sur¬ 
rounded  by  boundary  layers,  it  is  necessary  to  accurately 
approximate  the  inertia  terms  only  in  these  thin  regions. 

The  "modified  Oseen  technique"  replaces  the  convective 

velocity  there  by  a  given  suitable  average  velocity  and  there- 

1 8 

by  linearizes  the  Navier-Stokes  equations.  Weinbaum  has 
shown  (in  a  cylindrical  geometry)  that  this  approximation 
leads  to  a  smooth  transition  to  the  "Stokes"  solution  with  de¬ 
creasing  Reynolds  number. 


II.  Free  Shear  Layer  Approximation  and  Matching  Conditions 

On  the  DSL,  for  m  =  p  =  M  =  I,  the  matching  is  carried  out  at  one 
point,  chosen  to  be  r  =  R/2  (where  ~  prescribed  that 


(i) 

(ii) 

(iii) 


["dsl]*  =  0 
dsl]_  ~ 
Pdsl]* 


=  0 
2 


at  r  =  R/2 
at  r  =  R/  2 
at  r  =  R/2 


(21) 


Assuming  [l  -  Up+, 


we  require  that  Uj^  =  Uj^  . 

The  pressure  continuity  can  only  be  satisfied  approximately  since 
the  dividing  streamline  shape  is  given,  and  only  (i)  and  (ii)  will  be  satisfied 
here.  Then  u^^  is  determined  by  [(21)-(ii)]; 
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(22) 


Using  the  observation  made  before  concerning  free  shear  layer  solutions,  it 
is  assumed  that  the  velocity  profile  between  DSL  and  edge  of  the  free  shear 
layer  is  linear  and  the  shear  is  given  by 


D 


h  “  "d] 


From  Eq.  (20)  we  may  compute  T  : 

_  4ab 

Td  =  cosh  ap  sin  bp  +  sinh  ap  cos  bp] 

From  Eq.  (22)  and  Fig.  1, 


1  +  N 


,  S  I  Ism  a 


where 


[z  cosh  ap  sin  bp  +  sinh  ap  cos  bp], 


We  will  assume  6^  to  be  the  thickness  of  the  expanded  boundary 
layer  at  the  shoulder,  thus  neglecting  the  shear  layer  calculations.  This 
approximation  is  certainly  valid  for  thick  shoulder  boundary  layer  thick¬ 
nesses,  and  the  important  interaction  is  the  effect  of  Uj^  on  . 

III.  Shoulder  Expansion 

Since  an  exact  solution  of  the  corner  expansion  of  a  supersonic 
boimdary  layer  is  not  available,  the  relatively  simple  streamtube  method 
will  be  employed  to  estimate  ^2'^^z'  -^ssmning  an  isentropic  expansion  of 
each  streamtube  to  pressure  P2(determined  by  turning  the  inviscid  flow 
through  an  angle  v  =  a  +  p),  and  conserving  mass  in  streamtubes,  it  is 
possible  to  calculate  the  velocities  and  streamline  locations  after  expansion 
(details  of  this  calculation  are  given  in  Appendix  I).  Specific  results  are 
obtained  with  the  good  approximation  for  hypersonic  boundary  layers  of  a 
linear  velocity  profile  and  Busemann  Integral  enthalpy  profile.  The  boundary 
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1  3 

layer  thickness  before  expansion  is  taken  from  Truitt  for  insulated  walls, 

14 

and  it  is  seen  in  Schlichting  (for  representative  Mach  numbers)  that  cold 
wall  thicknesses  are  well  within  a  factor  of  two  of  the  adiabatic  results.  * 


IV.  Recompression  Process 

The  final  matching  condition  consists  of  equating  the  dividing  stream¬ 
line  pressure  at  the  stagnation  point  to  the  local  value  of  the  inviscid  pres¬ 
sure  field.  Since  the  inviscid  flow  is  not,  in  general,  horizontal  at  this 
location,  the  pressure  is  less  than  the  neck  pressure.  Assuming  isentropic 
processes  exterior  to  the  free  shear  layer, 


H 


(26) 


Then,  Pgj.^g  =  where  q  depends  on  the  "re -attachment"  process  in 

the  neck.  A  typical  value  calculated  by  Reeves^^gives  q  .  6.  It  is  now 
assumed  that  the  total  pressure  on  the  dividing  streamline  is  approximately 
constant  during  re  compression,  and  is  therefore  equal  to  at  the  matching 
point  (r  =  1/2).  This  assumption,  discussed  in  a  previous  section,  physic¬ 
ally  consists  of  a  division  of  the  dividing  streamline  into  two  parts:  a  region 
of  acceleration  by  shear  forces,  in  which  the  free  shear  layer  does  work  to 
increase  the  kinetic  energy  of  the  DSL,  and  a  region  of  isentropic  recom¬ 
pression,  in  which  further  pressure  rise  due  to  shear  forces  is  relatively 
insignificant. 

We  may  therefore  write 


1  + 


-  1  K/f  2 

2  ^3 


y/y- 


=  r)P 


H 


(27) 


It  is  to  be  noted  that  q  is  dependent  on  the  length  of  the  recompression  re¬ 
gion  in  the  free  shear  layer,  and  can  be  expected  to  decrease  as  Reynolds 
number  decreases. 


V.  Solution  of  the  Energy  Equation  in  the  Recirculation  Region 

Consistent  with  approximations  made  in  the  solution  of  the  momentum 
equations,  we  will  assume  the  density  to  be  approximately  constant  in  the 
low  Mach  number  recirculation  region.  We  may  then  neglect  friction  and 
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compression  work  in  the  energy  equation  and  write: 


where 


P6V*Vh  =  V  h 
V  =  iu  +  jv 

h  =  h^hj,  and  Pe-^  PrRe 


The  boundary  conditions  are: 


(i) 

h 

=  ffF) 

on 

0=  P 

(ii) 

3h 

Te 

=  0 

on 

0=0 

(iii) 

h 

=  f(l) 

on 

r  =  1 

where  f(r)  is  to  be  specified.  To  obtain  homogeneous  boundary  conditions, 

8Ct 

h  =  f(7)+h(r.  0)  (30) 


and  obtain 

-P^  u||  +  ^  11+  V^h  =  G(P<r,  0)  (31) 

.  £i 

where  G  =  +Pe  u  f  —  —  —  f”  is  a  known  function. 

r  _» 

The  resulting  nonhomogeneous  equation  is  linear,  since  V  has  been 

obtained.  An  expansion  of  h  in  Pe’-niimber  is  the  simplest  approach  to  a 

solution,  but  cannot  be  expected  to  converge  for  Reynolds  numbers  as  high 

as  those  found  valid  for  the  momentum  equation.  This  is  because  isotherms 

are  generally  expected  to  be  non -parallel  to  streamlines  due  to  the  wall  and 

axis  boundary  conditions,  and  PeV*  Vh  is  non -negligible  for  flight  Pe- 

numbers  (the  equivalent  term  in  the  vorticity  diffusion  equation  was 

Re  'V'VCiJ,  which  can  be  small  due  to  the  nearer  coincidence  of  vorticity 

9 

contours  and  streamlines  ).  We  must  therefore  retain  the  convective  terms 

and  deal  with  the  variable  coefficients  in  some  approximate  manner.  The 

procedure  adopted  is  to  solve  the  equivalent  "Galerkin"  formulation  of  the 

12 

problem  (see  Kantorovich  ). 
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Let 


oo  00 


(32) 


with  the  taken  to  be  a  complete  set  of  functions  satisfying  the  boundary 

conditions: 

(i)  ■h(0,  6)  =  0 

(ii)  h(l.e)  =  0 


(33) 


The  functions  =  sin  j  tt  r  are  therefore  suitable.  The  boundarv^  conditi 
(iii)  h(r,  (3)  =  0 


ions 


(iv)  ||(Y,0)  =  0 

are  satisfied  by  even  functions  comprising  a  "relatively  complete"  set; 


(34) 
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(35) 


gk(0)  =  (1  - 

Thus,  the  N/M  approximation  is  given  by 

M  N 

=  (l-eVp^)  22  22  a  e^^’^-^^sinjirr  (36) 

j  =  1  k  =  I  •’ 

The  solution  (the  a^j^' s)will  be  determined  from  the  equivalent  Galerkin 
formulation: 


Jq  ^  , 


gj^rdrd0  =  0 

for  j  =  1,  2 .  M 

k  =  1,  2,  3,  .  .  .  N 


(37) 


with 


u«»  =  16t2(1_„2 


v*°'  =  _  16r^(l  -r)(3 -5t) 

F(0)  =  z  sinh  a0  cos  b0  —  cosh  a0  sin  hQ 


from  Eqs.  (2)  and  (20) 
with  u  II  r^ 


and 
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/Pb^3^  \  / sin  a\ 

\  J  V^in  pj  U3 


We  have  assumed  the  Prandtl  number  to  be  constant  and  pg  to  be  approxi¬ 
mately  Pj/RT^  in  this  first  iteration.  To  simplify  the  integration  here, 
and  to  provide  polynomial  coefficients  for  a  possible  series  solution,  F(0) 
and  F'(0)  will  be  approximated  as  follows: 


F(0)  S  m0(l  -  0/p) 


F'(0)  S 


s  r  - 

[1  . 


1  =  rja^d -1 


It  is  seen  in  Figs.  (2)  and  (3)  that  F'(0)  is  very  well  represented  with 
«.  55  over  a  large  range  of  p.  The  approximation  of  F(6)  should  be  most 
accurate  near  0  =  p,  and  can  be  poorly  represented  near  0=0,  where 
3h/9  0-*O.  Thus,  m  is  required  to  minimize  the  average  weighted  error 
in  the  interval  (o,  p): 

I  e  |m0(l  -  0/p)  -F{0)j  d0  =  0  (40) 

Jr\ 


a  typical  solution  is  given  in  Fig.  (4),  and  m  is  given  in  Appendix  II. 

We  expect  that  thermal  boundary  layers  will  form  at  even  moderate 
P^-numbers,  and  the  approximate  solution  must  reflect  this  behavior.  Con 
vection  is  expected  to  decrease  temperatures  on  the  dividing  streamlines 
and  increase  them  on  the  axis,  indicating  that  a  reasonable  approximation 
may  be  obtained  with  M  =  1  and  N  =  2.  Proceeding  from  Eq.  (36) 


T  (2)  _ 

hi 


0Vp^)[aii 


and  we  must  solve  the  equations: 

^  j  [lCRj<^^)  -  g]  (1  -  0Vp 

^  -  g]  (1  -  aVp 


')simrrrdrd0  =  0 


0Vp^)0^  sin7rTTdrd0  =  0 
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v/here  = 

,  a' 


.  a%<^'  .r,„,a-./^'^,(o,  a^/a)| 

-1^  i-,r  ¥Tr  [  — 

2  2  r  2  2  2  —1 

=  —  (1  —  0  /p  )  a^  jff  sinTrF  +  a-j2®  sin^r 

1  2  2  r  —  2  T 

+  3  (1  —  0  /p  )  |a^  j^7rcos;rr  +  ^^^2^  cos:r7j 

1  r  ^^11  .  _  .  -  120^x1 

r  L  P  \  P  /J 


+  1  L  sin.-  +  a,  ,sin.T  (z  -  ^^-^1 

-Pej  16r^(l  -  r)^  [a^(l  -  0^/p^)  “  l]  [^  “  ^^/P^]  [a^^fl-cosTrr  +  Trai^^e^cosn^ 
-I6^(l-r){3-5r)  p.0{l_^p)  |1 sinTrr  +  Ze3^^2^imrr(^l 

To  proceed  further,  we  must  specify  f{T).  This  is  done  by  assuming  a 
linear  change  from  h^  to  hj^,  followed  by  a  constant  total  enthalpy  recom¬ 
pression.  These  assumptions  are  consistent  with  those  employed  in  solu¬ 
tion  of  the  momentum  equations  (Pj(o,  P)  «  Pj(l/2,  6)): 


for  .  5  <  r  <  1.0 


!  1  +  2^-  (ujj^  —  u^^)  for  0  <  r  <  .  5 


where  u*  =  u„u^°\  the  dimensional  velocity;  thus 

5|6  \J 


■h 

J  -  1  for  .  5  < 

"d 

L  ,  J  ,  u 


r  <  1. 0 
2 


f"(f)  = 


I- 512r^(l  -f)  (1  -2f) -5^  :  0<r<.5 

‘  D 

0  for  .  5  <  r  <  1.0 
!  ^2 

-  512r^(l  -  r)^(3  -  Ur  +  14r^) 

D 


The  Galerkin  equations  then  reduce  to: 
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ajjAj  +  Rj 


®11®1  +  ^2^2  =  ^2 


where  Aj,  K^,  B2,  Rj  and  are  given  in  Appendix  II. 

Solving  for  a^j  and  a^^i 


~7sjr 


•*2  “ 


®2-®i  jr 


and  it  is  noted  that  and  R2  depend  on  h^/hj^  and  Uj^/hj^, 
Now,  in  dimensional  terms. 


1  *'"♦1  "D  1  1 

7  "Sd"  ""  TT  7  ~?F" 

which,  at  F  =  1/2  and  0  =  p, 

4h„  , 


VI.  Free  Shear  Layer  Solution  of  the  Energy  Equation 


We  assume  that  the  enthalpy  profile  is  "locally  similar"  and  can  be 

14 

approximated  in  the  matching  plane  by  the  Busemann  Integral  solution 
(assuming  a  constant  pressure  free  shear  layer  with  Pr  =  1.0):  i.  e.  , 

H  =  Au  +  B,  with  the  boundary  conditions: 

(i)  H(u-)  =  H- 

^  ^  /  r>  /\  % 


Thus, 


(i) 

H(u3) 

(ii) 

H(Uj5) 

A  = 

Hj-Hd 

U3  -Ujj 
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and 


ah 


9u 


{A  -u) 


using  the  profile  of  Eq.  (23) 

«3 


ah 


(51) 


where  the  approximation  n  ^/2  has  been  used  and  the  definitions  made: 


(52) 


VII. 


Matching  Conditions 

The  matching  conditions  at  r  =  1/2,  6  --  |3  are  simply: 


(i)  [hp]^  =  0 


(53) 


Thus,  the  enthalpy  gradient  is  continuous  and  we  obtain: 


To  solve  for  hp,  we  must  write  I^  and  I^(integrals  which  appear  in 
Eq,  (42);  see  Appendix  II)  in  the  form: 
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and  define  h 


Now  defining 


32Q,Pe'(a^l,-Ij)-£l3  =  Nj 


1.64  X  10^Q2(a^Ig -Ig) -.46lg  = 
1.64X  lO^Q2(a^I, 46I„  = 


we  may  solve  for  hj^ 


.“‘,,2  ,  4(«3\ 

+FVir; - '- - '6-“--b;(a\/s7)— ^ - 


l(^) 


-2p  B 


. . 


We  may  now  calculate  u,  (subscripts  now  refer  to  the  iteration -number) 

from  Eq.  (25)  for  a  trial  p.  Inserting  u,  into  the  above,  '  '  is  known, 

#  1 
and  u_  may  be  calculated  from  Eq.  (27).  A  new  trial  p  is  selected  until 
♦  T  .  o  *  * 

u^^fUj,  starting  from  P  =  1  (u^  decreases  with  p  and  u^  increases  with  P; 

a  convergent  iteration  results). 

VIII.  Busemann  Integral  Solution  in  the  Recirculation  Region 

A  still  simpler  approximate  solution  may  be  obtained  by  extending 
the  free  shear  layer  solution  into  the  recirculation  vortex,  taking  as  a  new 
boundary  condition; 

H  =  h  at  u  =  0 
core 

hcore  best  approximated  by  assuming 

h  w  h  ,,  (56) 

core  wall  '  ' 
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(this  assumption  implies  the  absence  of  a  thermal  boundary  layer  on  the 
base  wall). 

Then,  T_  T  /  t\.  /  ,\  . 

Tf—  =  +  (l  “  T^"j  “  ^  u{l  -u)  (5 


which  replaces  Eq.  (55).  Actually,  the  solution  to  Eqs.  (25),  (27)  and  (57) 
may  be  solved  directly  by  assuming  (3,  solving  Eq.  (27),  with  Eq.  (57),  for 
u(a  quadratic  equation),  and  Eq.  (25)  for  Re^g.  For  the  sake  of  uniformity 
of  the  program,  however,  Tj^/T^  was  used  in  the  iteration  described.  The 
results  of  this  calculation  will  be  discussed  together  with  those  of  the  pre¬ 
vious  section. 

To  relate  the  remaining  unknowns  to  flight  conditions,  it  is  necessary 
to  write  down  the  inviscid  flow  field  solutions.  Given  the  velocity  and  al¬ 
titude,  ambient  conditions  are  known  and  post-shock  conditions  are  computed 
as  follows; 

piU.S 

M|=u,/aj,  Rej3=-_  (58) 


=  »,A| 


Shock  angle,  6,  is  obtained  from 
M,  ^sin^S  —  I  =  (■ 


-  1  =  (iji)  M, 


Then, 

and  with  M 


2  sinfi  sing 
cos(6  —  a) 


2  •  2. 
sin  0 


-1  1  +  — Mj  sin  6 

sin^(6  —  a)  sin^fi 


j  =  M.sinS,  standard  tables  are  used  to  compute  p,,p-5,  T.,,  a_ 
n  *  -tf  2  ^2  2  2 


and  U2.  We  assume  the  viscosity  law  Pg/M]  =  •  The  Prandtl- 

Meyer  Function  determines  M„  (required  for  calculation  of  P„)  and  M, 
from  an  expansion  of  a  and  a  +  P  degrees,  respectively.  Isentropic  equations 

give  , 


iy/y-i 


T 

3 

1  + 

Z  c 

T 

2 

1  + 

2  3 

P, 

*2 

TFJ-  = 
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From  the  integral  solution  mentioned  previously, 
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«2 

■sr 


where 


Re 


2S 


T575T^  ^®2s 

Pz 


1/2 


(62) 


(63) 


JL 

m. 


VI  +  m^ 

- —In  (n/  1  +  m_  +  'v/irTT  ) 

"*2  =  (^i^)  “2^ 

The  characteristic  Reynolds  number  for  the  recirculation  regions  is 


Re 


B 


w 


where 


(64) 


Of  course,  T^,  o  and  S  are  given  data  and  together  with  u^  and  altitude  con¬ 
stitute  the  independent  variables  in  this  analysis. 


IX.  Numerical  Results 

Calculations  have  been  carried  out  for  the  conditions  of  the  wind- 

tunnel  wedge  experiment  of  Dewey^^;  Mj  =  6,  a  =  15°,  S  =  .  0483  ft  and 

=  300°F.  The  base  pressure  was  calculated  over  the  Reynolds  number 

range  of  the  experiment  for  various  values  of  t|,  the  only  free  parameter 

in  the  solution.  The  theoretical  base  pressure  decreases  with  t^.  due  to  the 

reduction  of  the  dividing  streamline  stagnation  pressure  (Eq.  (27)).  Two 

results  are  shown  in  Fig.  5,  from  which  it  is  apparent  that  a  constant  value 

of  q  yields  a  reasonable  approximation  to  the  data.  The  value  q  =  0.  6  was 
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selected  based  on  the  theoretical  work  of  Reeves  (for  Re_„  3720). 

NlifOK 

It  is  seen,  however,  that  q  =  .4  results  in  a  theoretical  curve  that  is  never 
more  than  20%  from  the  data.  This  would  indicate  that  the  recompression 


process  is  rather  gradual,  with  the  neck  displaced  significantly  downstream 
of  the  DSL  stagnation  point.  Additional  calculations  were  carried  out  using  a, 
Busemann  Integral  approximation,  Eq.  (57),  It  is  seen  that  some  variation 
from  the  more  exact  solution  occurs,  due  to  the  core  temperature  being 
held  fixed  at  the  wall  temperature.  To  more  precisely  fix  the  stagnation 
enthalpy  and  to  carefully  investigate  the  effect  of  variable  wall  temperature, 
the  Gale rkin  solution  to  the  recirculation  region  enthalpy  field  is  a  necessary 
complication. 

A  series  of  calculations  was  performed  to  provide  a  rough  under¬ 
standing  of  the  full-scale  behavior  of  the  near  wake.  Four  representative 
wedges  are  considered:  S  =  1  ft  and  10  ft  and  h^  =  .  05,  .  15  (typical  of  low 
temperature  and  high  temperature  ablators).  The  wake-angle,  base  pres¬ 
sure  and  stagnation  enthalpy  were  calculated  at  various  free -stream  velo¬ 
cities  and  altitudes.  These  results  are  plotted  in  Figs.  6-8.  In  general, 
stagnation  enthalpies  decrease  significantly  with  increasing  altitude.  The 

3|t 

DSL  velocity  ratio,  u,  always  decreases  with  increasing  altitude,  while  hj^ 
slowly  increases  and  then  decreases  (at  approximately  200,  000  ft  for 
Uj  =  20,  000  ft/sec).  The  latter  effect  is  .due  to  the  competition  between 
thermal  conduction  and  convection  and  the  variation  of  the  stagnation  en¬ 
thalpy:  at  the  lower  altitudes,  the  stagnation  region  is  hot,  but  convection 
keeps  the  dividing  streamline  enthalpy  down.  As  the  altitude  is  increased, 
conduction  begins  to  become  dominant  and  increases  h^^slightly.  However, 
at  still  higher  altitudes,  the  conversion  of  kinetic  to  thermal  energy  during 
recompression  is  markedly  reduced  and  enthalpies  everywhere  begin  to 
drop. 

The  effect  of  increasing  wall  temperature  is  qualitatively  as  expected. 
Significantly,  it  is  seen  in  Fig.  8  that  stagnation  enthalpies  of  "cold"  ten -foot 
bodies  and  "hot"  one -foot  bodies  are  comparable.  At  equal  wall  tempera¬ 
tures,  however,  the  larger  body  has  a  stagnation  enthalpy  that  is  typically 
25%  higher.  Furthermore,  the  calculated  base  pressures  are  50%  higher 

with  a  relatively  hot  wall,  which  qualitatively  checks  the  experimental  ob- 

21 

servations  of  Kurzweg. 

From  Fig.  6,  we  observe  that  the  near  wake  "dimensions"  scale 
roughly  as  the  body  size  for  low  and  intermediate  altitudes.  At  the  higher 
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altitudes,  a  small  body  exhibits  a  proportionately  larger  base  region,  but 
this  never  approaches  the  size  of  the  near  wake  of  the  body  used  for  com¬ 
parison  (additional  calculations  show  that  the  wake  angle  increases  slowly 
with  increasing  wedge  angle). 

The  effect  of  free -stream  velocity  is  exhibited  in  Fig.  9.  Apparently, 
beyond  10,  000  ft/ sec  there  is  very  little  influence  of  this  parameter  on  the 
stagnation  enthalpy  ratio.  At  velocities  below  this,  the  greatest  effects 
occur  at  high  altitudes,  indicating  the  small  influence  of  variable  velocity 
during  a  ballistic  trajectory  (assuming  these  results  to  be  qualitatively 
correct  for  a  three-dimensional  body).  Similar  remarks  apply  to  the  wake 
angle  and  base  pressure  plotted  in  Figs.  6  and  7, 

The  effect  of  changing  the  only  free  parameter,  r|,  is  demonstrated 
to  be  small  in  Fig.  10.  Because  of  the  improved  agreement  with  the  experi¬ 
mental  base  pressures  at  q  =  .  4,  this  value  was  assumed  in  the  above  cal¬ 
culations. 

(2) 

It  was  also  found  that  the  approximation  solved,  h^  ,  was  suffi¬ 
ciently  general  to  allow  a  boundary  layer  to  develop  along  the  dividing 
streamline.  Enthalpy  profiles  at  two  altitudes  are  shown  in  Fig.  11,  where 
the  perturbation  from  the  boundary  value  is  seen  to  be  only  a  few  percent. 

It  would  be  possible,  of  course,  to  accurately  calculate  the  wall  boundary 
layer  with  more  ttrms  in  the  radial  direction.  From  the  present  computa¬ 
tion,  only  a  rough  estimate  of  the  base  heat  transfer  can  be  obtained; 

9h  (2)  ^ 

— ^ -0^/(3^)  (ajj  +  aj20^)7rcos7rr] 


and,  since 


^  b  = 


q\  (^>(0)  = 

“l 


F=  1 


,  the  stagnation  rate  is 


r 


where  Q  0  implies  heat  transfer  to  the  base.  The  variation  of  Q  (0)/ 
with  altitude  is  shown  in  Fig.  12,  where  it  is  always  seen  to  be  positive  for 
cold  wall  conditions  (h^  =  .  05)  and  negative  for  a  relatively  hot  wall  (h^=  .  15). 
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Fig.  9  Stagnation  Enthalpy  vs  Free  Stream  Velocity  at  Various 

Altitudes 


Fig.  10  Effect  of  t] -parameter  on  Stagnation  Enthalpy  at  Various 

Altitude  s 
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A2639 


Fig.  11 


Fig.  12 


e//J 


Static  Enthalpy  Profiles  at  r  =  .  5,  Normalized  to 


Heat  Transfer  Rate  at  Base  Wall  Stagnation  Point, 
Normalized  vs  Altitude 
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X.  Summary  and  Extensions 


A  tractable  and  physically  correct  model  of  the  two-dimensional 

near  wake  has  been  analyzed.  The  flow  field  was  divided  into  subregions, 

including  the  recirculation  regions,  which  were  matched  along  mutual  bound- 
th 

aries.  The  zero  order  solution  of  the  velocity  field  was  obtained  from  a 
series  expansion  of  the  streamf unction  in  local  Reynolds  number.  Calcula¬ 
tion  of  the  uncoupled  enthalpy  field  was  carried  out  with  a  Galerkin  procedure, 
with  these  results  valid  for  arbitrary  Reynolds  number.  Numerical  results 
were  obtained  for  the  conditions  of  a  wind-tunnel  wedge  experiment,  and 
agreement  of  the  base  pressure  variation  with  Reynolds  number  was  satis - 
factory.  Previous  theories,  ’  it  should  be  noted,  have  yielded  results 
which  are  Reynolds  number -independent.  The  effect  of  body  size,  wall 
temperature  and  free -stream  velocity  and  density  on  the  dividing  stream¬ 
line  stagnation  enthalpy  was  also  investigated.  It  was  found  that  the  stag¬ 
nation  enthalpy  ratio  increased  with  decreasing  altitude  but  was  rather  insen¬ 
sitive  to  Mach  number  at  hypersonic  speeds.  An  increase  in  wall  tempera¬ 
ture  was  found  to  significantly  increase  the  stagnation  enthalpy  and  base 
pressure.  It  was  seen  from  calculations  for  two  bodies,  that  near  wake 
size  scales  directly  with  body  size.  Finally,  it  was  shown  that  the  stagna¬ 
tion  enthalpy  was  relatively  insensitive  to  changes  in  the  ratio  of  dividing 
streamline  stagnation  pressure  to  neck  pressure,  the  single  free  parameter 
in  the  solution.  While  simple  improvements  are  possible  in  the  two- 
dimensional  analysis  (solution  of  the  free  shear  layer,  multiple -point 
matching  with  the  recirculation  regions,  etc.  ),  it  is  anticipated  that  three- 
dimensional  effects  will  be  more  important.  The  concentration  of  stream¬ 
lines  in  the  stagnation  region,  growth  of  the  free  shear  layer  displacement 
thickness  (due  to  geometry  alone),  and  the  incompletely  understood  recom¬ 
pression  process  external  to  the  dividing  streamline  are  expected  to  make 
quantitative  changes  in  the  near  wake  solution.  While  the  results  presented 
herein  are  expected  to  indicate  the  general  behavior  of  axisymmetric  near 
wake  properties,  it  must  be  emphasized  that  they  (and  all  other  presently 
available  theories)  are  strictly  applicable  to  two-dimensional  flow  fields. 
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APPENDIX  I. 


Expansion  of  the  Boundary  Layer 
The  geometry  of  this  problem  is  shown  below 


The  assumptions  are  as  follows: 

(i)  the  rapid  expansion  process  along  streamlines  is  essentially 
isentropic 

(ii)  each  streamline  reaches  the  pressure  P^  ,  as  determined  by 
turning  the  inviscid  flow  through  the  angle  v 

(iii)  the  initial  velocity  profile  is  linear  and  the  enthalpy  profile 
is  "similar”  =  Auj^  +  B) 

(iv)  Y  =  1.  4 

The  calculation  procedure  is  simply  to  expand  each  streamline  with  P^  and 
T^  held  constant,  and  calculate  the  velocity  and  density  after  expansion.  We 
then  conserve  mass  between  streamlines  to  locate  them  after  expansion,  and 
to  estimate  the  expanded  bovindary  layer  thickness.  Denoting  streamlines 
by  i  and  conditions  before  and  after  expansion  by  1  and  2,  respectively,  we 
can  write 


2. 


i  +  1 


n 


where  i  =  0,  1,  2,  .  .  .  .  n  and  y^  =  0 

o 
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Also, 


Pi. 

1 

_  U 

■ 

6 

1 

^1. 

yi. 

1 

1 

^  ■»? 

(equation  of  state) 


(given  profile) 
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M_  is  obtained  from  the  Prandtl -Meyer  Tables 

ForM,  >10,  we  may  use 
A  » 


It  should  first  be  pointed  out  that  under  the  given  boundary  conditions,  the 
expanded  flow  will  be  non-parallel.  We  must  therefore  assume  some  pro¬ 
cess  to  turn  the  inner  streamlines  back  to  the  direction  given  by  v.  A 
"lip  shock”  or  combination  of  expansions  and  shocks  is  probably  necessary, 
and  this  complicated  region  is  presently  under  study.  The  assumption  of 
constant  along  the  wall  is  also  an  obvious  inconsistency,  and  requires  a 

finite  velocity  to  exist  there  after  expansion.  However,  since  experiments 
.  19  .  . 

by  Hammit  indicate  the  reasonableness  of  this  approximation  away  from 

the  wall,  we  allow  "slip  velocities”  to  be  calculated  and  interpret  them  as 

22 

velocities  at  the  edge  of  a  high-shear  sub-layer.  A  typical  profile  is 
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shovm  in  Fig.  13  and  62/^1  plotted  against  all  the  independent  variables, 

M  i_  ,  T  /To  and  V  in  Figs.  14-15.  It  is  seen  that  a  significant  thickening 
can  occur,  and  the  results  are  certainly  limited  by  boundary  layer -shock 
wave  interaction.  The  most  rapid  changes  with  angle  occur  at  the  higher 
Mach  number  and  the  most  significant  Mach  number  effects  are  at  low  Mach 
number.  Increasing  wall  temperature  has  the  effect  of  decreasing  the 
thickening  and  at  a  temperature  ratio  of  0.  8,  the  Mj  =2  boundary  layer 
actually  becomes  thinner  for  small  angles  of  expansion  (this  is  actually  ob¬ 
served  in  Hammit's  Schlieren  photographs).  The  sub-program  detailed 
above  is  used  in  the  near  wake  solution  with  the  following  changes  in  nomen¬ 
clature 


a  +  p 
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Fig.  13 


i 

1 


i 


Prandtl -Meyer  Expansion  of  a  Boundary  Layer;  Velocity 


Profile  after  Expansion  at  Mj  =  2.00,  v  =  10  , 
=  0.  I  (V  =  l.  4) 

0 


Fig.  14  Prandtl -Meyer  Expansion  of  a  Boundary  Layer: 

6,/6,  vs.  V  and  M,  for  T  /T  =  0.  1 

W  Og 
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Fig.  15  Prandtl -Meyer  Expansion  of  a  Boundary  Layer: 

6,/6i  vs.  V  and  M,  for  T  /T  =  0.  9 
^1  ^  °5 
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^14  = 
^15  = 

^17 

^18 

^19 

a  = 
o 

Pe  = 

Pr  = 


13^^840 

l6p/35 

spV^is 

48pV3465 

34p^/5040 

zpV^s 

1.43 


l.O 


?  (t-M 


-“Tf 


[sin(bp)cosh(ap)(a^  —  +  2  ab) 

+  cos  (bp)sinh(ap)(a^  —  b^  —2  ab)  ] 
[sin(bp)sinh(ap){a  —  bz)  —  cos(bp)cosh(ap)(b  -!■  az)] 
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